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A study on the possibility of distinguishing new heavy Majorana neutrino models at LHC en- 
ergies is presented. The experimental confirmation of standard neutrinos with non-zero mass and 
the theoretical possibility of lepton number violation find a natural explanation when new heavy 
Majorana neutrinos exist. These new neutrinos appear in models with new right-handed singlets, in 
new doublets of some grand unified theories and left-right symmetrical models. It is expected that 
signals of new particles can be found at the CERN high-energy hadron collider (LHC). We present 
signatures and distributions that can indicate the theoretical origin of these new particles. The 
single and pair production of heavy Majorana neutrinos are calculated and the model dependence 
is discussed. Same-sign dileptons in the final state provide a clear signal for the Majorana nature of 
heavy neutrinos, since there is lepton number violation. Mass bounds on heavy Majorana neutrinos 
allowing model discrimination are estimated for three different LHC luminosities. 
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I. INTRODUCTION 

There are a lot of very convincing experimental data 
for neutrino oscillations and non-zero masses, The 
smallness of neutrino masses [1| is generally understood 
as a consequence of some see-saw mechanism. This brings 
the question of the possibility of new heavy neutrino 
states. So far none of these new states was experimentally 
observed [H, Q with masses up to Mjy ~ 100 GeV. For 
higher masses there are many suggestions of experimen- 
tal possibilities in the next high energy hadron-hadron 
colliders 0, in electron-positron linear accelerators 
[ESQ] an( l m ncutrinoless double-beta decay 0| . The 
properties of these new heavy states are the central point 
in many of the theoretical models proposed as extensions 
of the present Standard Model of elementary particle 
physics. Heavy Majorana neutrinos naturally generate 
B-L asymmetries and are good candidates for leptogene- 



sis 0, [l(J. Besides the masses and mixing angle values, 
an important point is the Majorana or Dirac nature of 
new heavy neutral states. This is directly connected with 
lepton number conservation and the general symmetries 
of any extended model. 

The authors already investigated some points in this 
direction for a single heavy Majorana production in 
hadron-hadron colliders [4[ . It was shown that this mech- 
anism could be more important than pair production and 
that like-sign dileptons can give a clear signature for this 
process. In the present paper we consider more gen- 
eral possibilities for a heavy neutral lepton production 
at hadronic colliders. If a new heavy neutral lepton is 
experimentally found, then a fundamental question to 
be answered will be its theoretical origin - singlet (VSM) 



Uj, doublet (VDM) [lj], mirror (FMFM) [13| and other 
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extensions. 

The main point of our work is to present a study of sig- 
natures and distributions that could, in principle, show 
the fundamental properties of possible new heavy neu- 
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the LHC Collaborations will concentrate theirs efforts on 
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new phenomena analyzing channels with low multiplic- 
ity and the simplest SM extensions. Since we presently 
have no signal for new additional gauge bosons, we will 
make the hypothesis that the dominant interactions of 
the new heavy neutrino states are given by the standard 
SUl(2) £g> Uy(1) gauge bosons. We can resume the new 
particle interactions in the neutral and charged current 
Lagrangians: 

Cnc = ~2^^ ll>1 ( gv ~ 9Al5)tpiZ^. (1) 

i 

and 



the mixing angles have been considered to be equal to 
the upper bound of the global approach. 

This study is organized as follows: in section II we 
present the total cross-sections for pair and single pro- 
duction for different models of heavy Majorana neutrinos 
at the proton-proton collisions at LHC energies, and the 
various distributions that were studied trying to distin- 
guish models, as well as the characteristics and number 
of events of the chosen processes. In section III, we dis- 
cuss the statistical treatment used to disentangle models 
and determine the theoretical origin of a heavy Majorana 
neutrino, if experimentally found. Finally, we present our 
conclusions. 
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where i,j stand for the usual fermions and the heavy 
Majorana neutrino (N). These models have only two pa- 
rameters: the new Majorana neutrino mass and a mixing 
angle, if we consider all the mixing angles with the same 
value. 

The general fermion mixing terms are given in [H] as 
well as the experimental bounds on the mixing angles 
(sin 9 m i x = s m ). For the purpose of the present investi- 
gation we resume these mixing angles in Table 1. The 
N — e — W vertex is of pure V — A type for the VSM, 
pure V + A type for the VDM and pure vector V or pure 
axial A for the FMFM. 

The study of the bounds on the mixing angles can 
be addressed in several ways. A general consequence of 
fermion mixing is that the standard model couplings of 
gauge bosons and fermions is to be multiplied by a fac- 
tor of the form cos 9 m ix ■ The universality of the differ- 
ent families coupling to the gauge bosons then imply a 
global upper bound on mixing angles that must all be 
very small, regardless of the studied model. The present 
value, with 95% C.L., for the upper bound on the mix- 
ing angle sin 2 8 m i X , is 0.0052, when considering them all 
equal. A detailed analysis of these calculations can be 
found in [4(. Another kind of analysis on mixing angles 
can be done for the experimental discovery potential at 
the LHC. In this case, besides the mixing angles one must 
add to the analysis at least one new unknown parameter 
- the heavy neutrino mass. This kind of approach was re- 
cently presented in reference FJ| for the V — A case. The 
possibility of new gauge interactions can add more un- 
known parameters to the analysis. The other V + A,V, A 
cases can be studied in a similar form. In this paper all 



II. TOTAL CROSS-SECTIONS AND 
DISTRIBUTIONS 

In our tree level calculations, we have implemented the 
models in the program CompHEP 16[ and have used 
the following cuts to take into account the detector ac- 
ceptance: \rj e \ < 2.5 (pseudo-rapidities of the final elec- 
trons); \t]w\ < 5 (pseudo-rapidity of the W); E e > 5 GeV 
(electron energies), and M ee > 5 GeV (invariant masses 
of the pairs of final electrons) . 

The total cross-section for single and pair production 
of heavy Majorana neutrinos at the LHC energy is given 
in Figure 1. As a consequence of the mixing terms, for 
the V—A case the pair production is strongly suppressed. 
Pair production is the dominant mechanism for pure vec- 
tor and axial couplings. We found no sensible difference 
between axial FMFM and vector FMFM for the total 
cross sections and for the distributions studied. So here- 
after we will mention FMFM for both cases. For the 
V ± A case, the dominant process is the single heavy 
neutrino production. Using a luminosity of 100 fb _1 , the 
upper bound on the Majorana mass is 500 GeV for the 
single production case. 

The next important point is the nature of the Ma- 
jorana couplings. Let us analyze first the real pro- 
duction of single Majorana neutrinos. For the process 
p+p — ► e+ + N + X the primary electron rapidity could 
indicate the V + A nature of the couplings, as shown in 
Figure 2. From these figures we can see clearly that there 
is no significant model dependence and the main differ- 
ences come from the total cross-section. 

Since the heavy Majorana neutrino can decay as N — > 
e T + VF ± , the most interesting, practical and clean signa- 
ture is the process p+p — ► e T + e T + W + X, with two 
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TABLE I: Mixing angles for W and Z couplings with new neutral leptons for VSM, VDM and FMFM. 




FIG. 1: Total cross section and number of events for p + p — > 
e~ + N + X andp+p-> N + N + X for VSM, VDM, FMFM 
at y/s = 14 TeV for C = 100 fb _1 , where TV stands for the 
new heavy Majorana neutrino. 

identical electrons or positrons in the final state. The ra- 
pidity of the final two electron momentum sum is given 
in Figure 3. In order to avoid double counting of the 
final identical particles we adopt to sum their momenta 
and calculate the rapidity of this sum. These are more 
realistic distributions, for practical purposes, but again 
there is no strong model dependence in their shapes. 

In Table 2 we show the expected number of events for 
the process p + p — ► e~ + e~ + W + + X at = 14 
TeV, for a luminosity of 100 fb , corresponding to one 
year of LHC operation at high luminosity regime. If we 
also consider the process p + p — ► e+ + e+ + W~ + X, 
the number of events in the table should be multiplied 
by two. 

Although not so visually distinctive , model separation 
can be achieved through a statistical treatment using the 
rapidity of the lowest px electron in the process p+p — > 
e T + e T + VF ± + X . The rapidity of highest px electron 
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FIG. 2: Rapidity distributions for the final electron momen- 
tum in p + p — > e~ + N + X for Mat = 100 GeV (top) and 
Mjv = 400 GeV (bottom) for VSM, VDM, FMFM at s/s = 14 
TeV. 

showed to be less sensitive to characterize the difference 
between the models. The lowest pr electron distributions 
obtained for M N = 100 GeV and M N = 200 GeV for C = 
300 fb -1 are showed in Figure 4 for the different models 
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p + p-> e + e " + W* + X 
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FIG. 3: Rapidity distributions for the sum of the final two 
electron momenta in p+p — > e~+e~+W + +X for Afjv = 200 
GeV (top) and M N = 400 GeV (bottom) for VSM, VDM, 
FMFM at s/s = 14 TeV. 
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FIG. 4: Rapidity distributions for the final e with lowest pr 
in p + p — > e~ + e" + W + + X for M N = 100 GeV (top) 
and M N = 200 GeV (bottom) with C = 300 fb" 1 for VSM, 
VDM, FMFM at y/s = 14 TeV. 



studied. One can see from Figure 4 (top) that it is not 
difficult to separate FMFM from the VSM and VDM, 
mainly due to the different cross sections, but the VSM 
and VDM separation is harder to do even for a Mjv = 100 
GeV. Figure 4 (bottom) shows that for M N = 200 GeV 
it is already impossible to distinguish among the models, 
unless one uses a statistical treatment. 
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TABLE II: Number of events expected for the process p + 
p — > e~ + e~ + W + + X at y/s = 14 TeV for VSM, 
VDM, and FMFM, for some Majorana neutrino masses and 
for C = 100 fb" 1 . 



III. MODEL DISTINCTION 

We performed a statistical treatment using an approx- 
imate x 2 -function for Poisson distributions, adapting the 
Eq.(25) of reference 15| in order to compare the his- 
tograms and to find the Majorana mass limits allowing 
model separation for three different luminosities. We use 
the following x 2 -function 



V" 2(rii - h) + (2m + 1) In 



/ 2rii 



1 



\2ki + l 



(3) 



where and hi are the contents of the i th bin of each of 
the two histograms. 

Our choice for Eq.(3) shows better results than the 
usual one [l| , even in the cases when the number of events 
is small, especially when there are some bins with zero 
contents. For large number of events, Eq.(3) behaves 
asymptotically like the least square method. It converges 
faster, with much smaller number of events, and with less 
oscillations to the true parameter values when used for 
fitting data as shown in 15j. In order to be more realis- 
tic it was included an efficiency of 80% for the electron 
identification and reconstruction. 

The events generated by CompHEP have been used 
to calculate the x 2 -f unc ti° n differences among the his- 
tograms for different models and masses. The calcula- 
tions have been done at the parton level and the final 
W disintegration and hadronization have not been per- 



formed, since we are using only information from the final 
electrons and they allow us a clean analysis. A clear sig- 
nature is to require two identical electrons and hadron 
jets in the final state. Due to the low number of events 
and to analyze the histograms fluctuations, we have re- 
peated the x 2 -function calculations several times using 
10 different sets of data for each model, mass and lumi- 
nosity and taken the mean values and their respective 
errors. This means that for each two model comparison, 
the x 2 -function calculation has been done 100 times in 
order to obtain the mean values and their respective er- 
rors. Since we have two identical electrons in the final 
state, we have chosen the rapidity of the electron with 
the lowest px to calculate the modified x 2 • This variable 
distribution was shown to reinforce the differences among 
the models under study. We have not worried about the 
backgrounds, since this process is forbidden in the Stan- 
dard Model and the background electrons coming from 
the decays of Standard Model particles have low energies. 

The results for the separation among the models can 
be seen in Figure 5. The straight horizontal solid line 
represents the value above which models can be discrim- 
inated with a confidence level larger than 95%. This 
value has been calculated for 50 degrees of freedom since 
we have used 50 bins in the < x 2 >-function calculations. 
In Figure 5 (top) we can see that for the VDM-FMFM 
combination, the models can be distinguished for Mn 
above 250 GeV, with data corresponding to one year of 
LHC operation at high luminosity. For a luminosity of 
300 fb -1 , VDM and FMFM can be separated if the heavy 
neutrino mass reaches 330 GeV. For 500 fb -1 , the sepa- 
ration between models can be achieved up to Mn ~ 400 
GeV. Similar mass bounds have been obtained for the 
separation between VSM and FMFM in Figure 5 (mid- 
dle). On the other hand, for the VSM- VDM separation, 
heavy Majorana mass limits were found not so high, as 
shown in Figure 5 (bottom). This difference between the 
models reflects the different couplings (V— A, V+A, pure 
V or pure A). The distinct total cross sections increase 
the differences in the rapidity distributions. For FMFM, 
no significant distinction was found between pure V and 
pure A couplings, as said before. In our calculations, 
FMFM was considered as pure V couplings. 

In order to better quantify the above results, p- values 
are also calculated. The p-value, in this paper, is the 
probability of a hypothesis given by a histogram Ha be 
incompatible with another hypothesis given by histogram 
Hb ■ If the two histograms are strongly incompatible then 
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0, in the other hand if they are strongly compatible 

1. The p- values are calculated according to 



/>oo 

P = / f(z,n d )dz 
J <x 2 > 



(4) 



where f(z, rid) is the x 2 probability distribution function 
and rid is the number of degrees of freedom, in this study 
rid = 50. < x 2 > are the values obtained when compar- 
ing the histograms and shown in Figures 5. The Tables 
III-V show the p- values obtained for different masses and 
luminosities. If we consider a model separation with 95% 
C.L. then p < 0.05. If the data at LHC show a compat- 
ibility with one of the three models described here, then 
it will be possible to distinguish the other two according 
to the Tables III-V. 
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FIG. 5: Separation between the models VDM-FMFM (top), 
VSM-FMFM (middle) and VSM- VDM (bottom), the hori- 
zontal line of each figure corresponds to 95% C.L. based on 
the < x 2 > with 50 degrees of freedom. Note the logarithmic 
vertical scale for the top two figures. 
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M N 
{GeV) 


c = 100/&- 1 


C = 300/6" 1 


C = 500/6" 1 


100 


0.102651 ■ lO" 120 


0.275233 ■ 10~ 387 


0.42576 ■ lO" 660 


200 


0.588073 • 10~ 4 


0.691069 • 10~ 9 


0.29823 ■ 10~ 15 


250 


0.386516 • lO" 1 


0.664087 ■ 10~ 6 


0.95806 ■ 10~ 9 


300 


0.931739 


0.360976 ■ 10^2 


0.891443 ■ 10~ 6 


400 


0.999999 


0.913923 


0.157242 



TABLE III: p-values for different masses and luminosities 
when comparing VDM and FMFM 



Mjv 
{GeV) 


C = 100/6- 1 


C = 300/6- 1 


C = 500/6- 1 


100 


0.304084 • 10" 203 


0.233734 • 10~ 633 


0.434935 • 10" 1073 


200 


0.224681 ■ lO" 11 


0.108177 ■ 10~ 47 


0.837551 • 10~ 81 


250 


0.289768 • 10~ 2 


0.247785 ■ 10~ 20 


0.231413 • 10~ 36 


300 


0.843227 


0.119454 ■ 10~ 5 


0.266346 • 10~ 14 


400 


0.999999 


0.782155 


0.270090 ■ 10- 1 



TABLE IV: Same as Table III but for VSM and FMFM 



Mjv 
(GeV) 


C = 100/6- 1 


C = 300/6" 1 


C = 500/6" 1 


100 


0.958060 • 10~ 9 


0.105196 ■ 10~ 22 


0.747901 ■ lO" 40 


200 


0.728142 • 10~ 3 


0.194823 ■ 10~ 12 


0.7882270 ■ lO" 20 


250 
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0.664087 ■ 10~ fi 


0.184826 ■ 10~ 9 
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0.998807 


0.157242 
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0.999999 


0.994980 


0.6345808 



TABLE V: Same as Table III but for VSM and VDM 
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IV. CONCLUSIONS 

This work suggests that the process p + p — ► e T + 
is an adequate, relatively simple and clear 
channel to look for a heavy Majorana neutrino and also 
to study its origin at hadronic colliders, since it violates 
leptonic number conservation rule and has, a priori, no 
background from the Standard Model. For an annual lu- 
minosity of 100 fb _1 this channel allows an investigation 
of the existence of a Majorana neutrino up to a mass of 
500 GeV. Furthermore, we showed that, once experimen- 
tally found, the theoretical origin of the Majorana neu- 
trino can be determined using the rapidity distribution 
of the final electron with the lowest pt, provided there 
are enough data. Although the differences among the 
rapidity distributions for the several models studied are 
not so great, we showed that a model separation could 
be achieved using a special statistical treatment, until 
a mass of the order of 300 - 350 GeV for 300 ib" 1 of 
data, at the parton level. The use of the lowest pr elec- 



tron rapidity distribution enhances the difference among 
the models, total cross sections and thus allows model 
discrimination until higher Majorana masses. A treat- 
ment, considering the hadronic jets resulting from the W 
disintegration, could lead to model separation for higher 
heavy neutrino Majorana masses. The LHC experimen- 
tal groups can do a more realistic treatment doing full de- 
tector simulation. Due to the few number of parameters 
involved in the models studied here and the clear non-SM 
channel, we believe this kind of experimental analysis can 
be done soon with data from the first " runs" . 
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